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ABSTRACT 

This paper presents a new study of the spectral index distribution of the supernova 
remnant (SNR) Puppis A. The nature of field compact sources is also investigated 
according to the measured spectral indices. This work is based on new observations of 
Puppis A and its surroundings performed with the Australia Telescope Compact Array 
in two configurations using the Compact Array Broad-band Backend centred at 1.75 
GHz. We find that the global spectral index of Puppis A is a = —0.563 ± 0.013. Local 
variations have been detected, however this global index represents well the bulk of the 
SNR. At the SE, we found a pattern of parallel strips with a flat spectrum compatible 
with small-scale filaments, although not correlated in detail. The easternmost filament 
agrees with the idea that the SNR shock front is interacting with an external cloud. 
There is no evidence of the previously suggested correlation between emissivity and 
spectral index. A number of compact features are proposed to be evolved clumps of 
ejecta based on their spectral indices, although dynamic measurements are needed to 
confirm this hypothesis. We estimate precise spectral indices for the five previously 
known field sources, two of which are found to be double (one of them, probably 
triple), and catalogue 40 new sources. In the light of these new determinations, the 
extragalactic nature previously accepted for some compact sources is now in doubt. 
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1 INTRODUCTION 

Puppis A is one of the most interesting supernova rem¬ 
nants (SNRs) in the Southern hemisphere. Its large extent 
(~ 50 arcmin in diameter) and high surface brightness pro¬ 
vides an unparalleled chance to study detailed structures. 
Together with Cassiopeia A and G292.0+1.8, it is one of 
the three oxygen-rich Galactic SNRs, where fast oxygen 
knots of still uncont aminated ejecta have been detected at 
optic al wavelengths JW inkler fc Kirshnen llOSSl l and in X- 
rays (iKatsuda et al.l[2008l l. Relative abundances of metal- 
rich ejecta measured on the NE portion are consistent with 
a progenitor with a mass in the range from 15 to 25 M© 
dHwang et al.ll20o'^ : IKatsuda et al.llioi^ 'l. 

The compact X-ray source RX J0822-4300, near the 
centre of the shell and hence dubbed ‘central compact ob- 
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ject’ (CCO), has been identified as the stell ar remnant 
left behind after the su pernova (SN) outburst dPetre et al.l 
19961:IZavlin et al]ll999l l. The explosion site was spotted by 


Winkler et al.l ( 198811 through measuring proper motions of 


ejecta knots and filaments. This positi on is fully consistent 
with the direction of the CCO velocity (iHui fc B eckeill2 006l: 
IWinkler fc Petrell2007l : iBecker et al.ll20ia r Castelletti et alJ 
( 2006tl propose that the CCO could be ejecting two opposite 
jets whose termination shocks appear in radio and X-rays as 
the eastern and western ‘ears’. 

The age of Puppis A is estimated to vary between 
3700 yr, based on th e proper motion of optical filaments 
d Winkler et al.l Il988ll . and 5200 yr, based on the veloc¬ 
ity of the CC O and its project ed distance to the ex¬ 
plosion centre d Becker et al .11201^ 1. The brightest, eastern 
rim of the shell is apparently inte r acting with a molec- 
ular cloud I Dubner fc Arnall 1 19881 : iRevnoso et al.l Il995l : 
IParon et al.l 1200^ 1. whose systemic velocity of ^ -|-14 km 
s“^ translates into a kinematic distance of 2.2 kpc. However, 
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IWoermann et all (l 200 (t) observed several pointings towards 
Puppis A and in the immediate vicinity in the four 18 cm 
lines of OH and concluded that the systemic velocity of the 
remnant is 7.6 km s“^ rather than +14 km s“^, which would 
bring the SNR about 1 kpc closer. 

In X-rays, Puppis A is one of the brightest SNRs, and 
the emission is mostly dominated by shocked interstella r 
medium (ISM: lHwang et aLll2005l . [200^ : lDubner et al.ll2013ll . 
Infrared (I R) observations car ried out with the Spitzer Space 
Telescope ijArendt et ahllioinl l show an extremely good cor¬ 
relation with X-ray emission, indicating that the thermal 
IR emission mostly arises from swept-up dust collisionally 
heated by the hot, shocked plasma. Puppis A has also 
been detected i n 7 -rays between 0 .2 and 100 GeV with the 


Fermi satellite llHewitt et al.l 20 

3), but not above 260 GeV 

(IH. E. S. S. Gollaboration et al.l 

2015|) . The lack of emission 


at TeV energies implies that a spectral break or cutoff would 
occur at 280 or 450 GeV, depending upon the power law 
model assumed. A break or radiative cutoff are expected if 
the acceleration of particles after the collision with an ex- 


ternal cloud has now ceased ||H. E. S. S. Gollaboration et al.l 

|2015l)- 

Puppis A was widely observed at radio-continuum 

wavelengths, both with single dish telescopes 

(Milne & Hilll 

Il969l: lKundulll97nl; lMilnelll97ll: iMilne et al.ll 

993h and in- 

terferometers dGreen 

I 1971 I: IMilne et al. Il983: 

Arendt et al.l 

I 1990 I: iDubner et al. 

1991: iGastelletti et al.l 

M061. Analv- 

sis of the polarized emission dKundu 

Il97r 

iMilnel 1 19721: 

iDickel & Milnelll976l: IMilne et al.lll99f 

) shows a predomi- 


nantly radial magnetic field, low polarization, and little de¬ 
polarization over most of the remnant. The first spectral 
index estimation of Puppis A was done by iMilne fc Hilll 
lll969ll who, based on the computed value of —0.5, con¬ 
cluded that this radio source was an SNR. Subsequent stud¬ 
ies refined the dete rmination of the bulk spectral index as 
Q = -0.52 + 0.03 ll^telletti et all l2006ln or a ~ -0.53 
l|Milnel Il97ll : IPubner et al.lll991^ . High-resolution studies 
have unveiled that the spectral index varies spatially ac ross 
the remnant llDubner et al.in~99ll : lGasteIletti et al.ll200^ '). In 
this paper, we revisit the spectral distribution of Puppis A 
applying a T-T plot method, which removes the problem of 
the zero level that can distort the results if the background 
emission is not properly accounted for. Our analysis is based 
on new wide band radio observations performed with the 
Compact Array Broad-band Backend (CABB) of the Aus¬ 
tralia Telescope Compact Array (ATCA). In addition, we 
re-analyse five previously known compact sources and list 
about 40 new compact sources in the field. 


2 OBSERVATIONS AND DATA REDUCTION 

Puppis A was observed simultaneously in the Hi 21 cm line, 
in four OH 18 cm lines, and in radio-continuum with the 
CABB of the ATCA in two 13 hr runs. The first run took 
place on 2012 May 20 using the array configuration EW 
352 (baselines varying from 31 to 352 m excluding the 6th 


^ Note that, as an in termediate step in the data processing, 
ICastelletti et all ll2006ll used two flux densities, at 86 and 408 
MHz, and found that a = —0.6. 


antenna), and the second, on 2012 July 6 in the 750A con¬ 
figuration (baselines from 77 to 735 m). Including the 6th 
antenna adds baselines between 3015 and 3750 m for the 
750A configuration, and between 4087 and 4439 for the EW 
352 configuration. The primary flux and bandpass calibrator 
was PKS 1934-638, while PKS 0823-500 was used for phase 
calibration. Due to the large angular size of Puppis A, ob¬ 
servations were made in mosaicking mode with 24 pointings 
in order to cover not only the SNR but also its surroundings. 
Following the Nyquist theorem, the pointings were separated 
by 19.6 arcmin to optimize the sampling at 1.4 GHz. The 
radio continuum was observed using a 2 GHz bandwidth di¬ 
vided into 2048 channels of 1 MHz width, centred at 1750 
MHz. The analysis of the Hi and OH lines will be developed 
in subsequent papers. 

The data re duction was carri ed out with the miriad 
software package llSault et al.ll 199^1 . The uv data were split 
in bands of 128 MHz. For each band, calibrator and source 
observations were inspected and flagged for interference and 
corrupted data. The standard data calibration routines were 
applied. We have not used the whole set of visibilities to 
construct a map because the presence of sources with very 
different spectral slopes within the observed field will intro¬ 
duce non-negligible errors due to bandwidth smearing in a 
frequency band as wide as 2 GHz. Therefore, we preferred 
to use at most three consecutive bands in a single image. 
Using a multifrequency technique, we constructed an image 
at 1.4 GHz merging the bands centred at 1302, 1430 and 
1558 MHz. The data were Fourier-transformed with uniform 
weighting to reduce sidelobes, and cleaned with MOSMEM. 
The resulting image, shown in Fig. [T] has a beam of 82.2 
X 50.6 arcsec, with a position angle of —0?55, and a noise 
level of 1.5 mjy beam“^, and is centred at 1395.6 MHz. A 
few residual sidelobe rings from the bright compact source 
at RA(2000)=8''20“35!6, Dec.(2000)=-43°0'23" could not 
be completely removed. The integrated flux density is com¬ 
puted to be 115 ± 8 Jy, where the error is dominated by the 
uncertainty in the background emission and the outermost 
contour level adopte d. The agreemen t betw een this value 
and that reported bv iGastelletti et al.l ll2006l ) indicates that 
at this frequency the ntscoverage, combined with the mo¬ 
saicking technique, samples all scale structures within Pup¬ 
pis A. 

To compute spectral indices, we applied a standard T-T 
plot method (see sni- This method requires that the two 
images to be compared, one at each frequency, be obtained 
from visibilities within a same uv range and convolved with 
the same beam. The first condition ensures that both im¬ 
ages sample features with similar scales. Adopting the first 
frequency to be 1.4 GHz and the second 2.5 GHz, visibilities 
were found to overlap in the range between 0.2 and 4 kA. 
We therefore constructed an additional image at 1.4 GHz 
following the same method as in Fig. [T]but using only this 
limited visibility range. The Mi)-filtered image thus obtained 
has a beam of 80.8 x 49.5 arcsec, with a position angle of 
—0?4. The image at 2.5 GHz was constructed merging the 
visibilities within the same uv range (0.2-4 kA) in the bands 
centred at 2326, 2454 and 2582 MHz. This image was con¬ 
volved with the same beam as the mia filtered image at 1.4 
GHz. To enhance the compact sources reducing the diffuse 
emission, two more sets of maps were formed with more re- 
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Figure 1. Image of Puppis A at 1.4 GHz. The beam size is 82.2 x 50.6 arcsec, with a position angle of -0?55 degrees. The flux density 
scale is shown at the right. The noise level is 1.5 mjy beam“^. The beam is plotted at the bottom left corner. The white numbers 1 to 
5 correspond to the compact sources discussed in E2] 


strictive uv ranges, from 0.5 to 4 kA and from 1 to 4 kA, 
respectively. 

Finally, to perform an accurate analysis of the compact 
sources, an image was constructed for each band applying a 
more restrictive visibility range, from 0.8 to 4 kA, to hlter out 
extended emission. These images were then used to compute 
flux densities and sizes, minimizing the uncertainties intro¬ 
duced by underlying background emission. The reduction in 
the number of visibilities due to the use of one single band 
and the restricted uv range results in a poor imaging and 
a higher noise, which makes the recognition of the weakest 
sources more difficult. Such sources were confirmed when 
they were detected in most frequency bands. 


3 RESULTS 

3.1 Spectral index computation 


To compute spectral indices minimizing the uncertainty in¬ 
troduced by the always arbitrary selection of the zero level, 
we have applied a T-T plot technique. A description of this 
method, widely used over sele c ted regions in extended ra¬ 
dio sources (e.g. Costain[ 1960l: Anderson_&_^Ajdrucld 

' e3l 


lUvamker et al ] |2004 |Leahvll2006l ?' iReynoso fc Green||2007l) 

can b e found in Turtle et al.l (1 19621 and Reich fc Reichl 
). Basically, the idea is to assume that the observed 


continuum emission at a given frequency T{i/) can be writ¬ 
ten as T{u) = Tslu) + To{i'), where Ts{u) is the source emis¬ 
sion and Tq{v) is an isotropic background component, which 



Figure 2. T-T plot between 1.4 and 2.5 GHz for the whole SNR 
Puppis A. Clip levels of 15 and 10 mJy beam“^, respectively, 
were applied. Background sources were removed. The fit yields a 
spectral index of a = —0.573 ± 0.004. 


includes an offset zero. Then, if Ts{v) oc with /3 constant, 
the brightness temperature at two different frequencies are 
linearly related as 

T{v,)^T{v2)['^y + C, (1) 

where C — Tq{vi) — Tq{v 2 ) {viIv 2 )^■ Thus, plotting T{vP) in 
terms of T{v 2 ) for all pixels within a given region, the value 
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Table 1. Catalogued compact sources 


Source 

number 

R,A (J2000) 
h m s 

Dec. (J2000) 

0 / // 

Spectral 

index^ 

Spectral 

index'^ 

Previous 

names 

1 

08 20 35.6 

-43 00 23 

-0.76±0.04 

-0.86±0.04 

G260.2-3.7H) 

2 

08 22 13.4 

-43 17 38 

-0.91±0.07 

-1.09±0.05 


3 

08 22 27.0 

-43 10 29 

-0.26±0.07 

-0.56±0.5 

ATPMN 1082226.8-431026^) 
2XMM 1082226.9-431026^) 

4 

08 22 45.7 

-43 12 38 

-0.75±0.15 

-0.94±0.08 

G260.7-3.7H) 

5 

08 25 30.7 

-42 46 57 

-0.92±0.04 

-1.09±0.09 

Cul 0823-426H) 


® T-T plot. ^ Slope of S(i^) versus u. 

Notes. I'll lPav et al.l lll972h : (21 iMcConnell et al.l ll2012h : fSf lWatson et all ll200Sll : (4) [sI3 |l973) 


of the brightness temperature spectral index /3 is readily ob¬ 
tained from the slope of the linear regression independently 
of the zero level of each of the two maps. 

In what follows, we will refer to the flux density spectral 
index a defined by 3(1/) oc i'°‘, which is related to the bright¬ 
ness temperature spectral index as a = /3 -|- 2. As the region 
where this method is applied becomes smaller, there is a bet¬ 
ter chance that both hypotheses, i.e. background isotropy 
and constant a, approach rea lity. A limitation of th i s tech - 
nique, already pointed out bv I Anderson fc RudnickI (Il993ll . 
is that spectral indices can only be measured at scales com¬ 
patible with the region analysed, while larger features are 
filtered out. 

The quoted uncertainties correspond to the for- 
m al statistical errors of the linear re gressions scaled by 
y/ (beam area/pixel area) (|Greerjri99flll . The scaling factor 
is due to the oversampling of the images, otherwise the sta¬ 
tistical errors in the linear fits may be underestimated since 
adjacent pixels are not independent a nd there are typically 
many pixels per beam (see for example iGaensler et al.|[l999l : 
iDohertv et al.l l2003h . The quoted uncertainties, therefore, 
account for oversampling and contain possible small-scale 
variations of the local background emission. In this estima¬ 
tion, the contribution of the flux density noise levels of each 
frequency map has not been taken into account. 

Although T-T plots are less accurate for extended 
sources, we used this method to estimate a global spectral 
index for Puppis A. We first plotted the nn-filtered image at 
2.5 GHz in terms of the same at 1.4 GHz (Fig. [2|). To avoid 
spurious contribution from emission unrelated to Puppis A, 
the brightest point sources projected within the SNR shell 
were removed in both images. The 1.4 and 2.5 GHz images 
were clipped at 15 and 10 mjy beam“^, respectively. A lin¬ 
ear fit yielded a spectral index of a = —0.573 ± 0.004, with 
a correlation coefficient of 0.9946. Inverting the frequencies, 
the result is a = —0.554 ± 0.004. Averaging both results, 
the global spectral index of Puppis A is a = —0.563 ±0.013, 
where the error encloses the two values obtained above. 

To map the spectral index distribution, we considered 
each pixel as the centre of a 7x 7 pixel box and fitted a T-T 
plot over the box between 1.4 and 2.5 GHz. The box was 
slid in two dimensions over the SNR at 1 pixel increments. 
We adopted clip levels of 20 and 15 mJy beam“^ respec¬ 
tively, and all linear fits where the correlation coefficients 
were lower than 0.950, or where the regression involved less 
than 6 pixels, were blanked out. This procedure was repeated 
for S'2.5 GHz against 5'i,4GHz and vice versa, and the average 


Table 2. Spectral variations for sources 3 and 4 


Flux density range 
at 1.4 GHz (mJy) 

Flux density range 
at 2.5 GHz (mJy) 

Spectral 

index 

5 ^ 45 

S 35 35 

-0.26±0.06 

3 45 ^ S ^ 100 

35 < S sS 80 

-0.5±0.3 

S ^ 150 

S ^ 120 

-0.14±0.20 

5 ^ 30 

S 35 20 

-0.7±0.2 

4 30 ^ 5 80 

20 S 50 

-1.0±0.4 

S ^ 80 

5 35 50 

-0.55±0.35 


between both maps was computed. The result is shown in 
Fig.m In good agreement with the global spectral index es¬ 
timated through the T-T plot shown in Fig. [5] the average 
spectral index of Puppis A obtained over Fig. Rafter removal 
of the background sources gives a = —0.59 ±0.22, where the 
larger error is due to spatial variations of the synchrotron 
spectra. 

3.2 Field compact sources 

A number of compact sources, either projected within the 
SNR shell or in the surroundings, are present in our maps. 
However, on ly five of them have been reported in previ¬ 
ous surveys (Day et ahl 1 19721 : ISled Il977l : iMilne et al.lll983l : 
iDubner et al. [l99lh . four of them superimposed on the SNR 
face. In Fig. [T] these sources have been labelle d with num¬ 
bers increasing with Right Ascension (note that iMilne et al.l 
1 19831. name them with letters A to E foll owin g the same or¬ 
dering ) . According to lMilne et al.l (Il983ll and lPubner et ahl 
11199111 . sources 1 to 4 have spectral indices close to —1.0, 
pointing to an extragalactic origin. Those indices are based 
on the direct comparison of integrated flux den sities at two 
frequencies. Source 5, catalogued as 0 823-4 26 bv ISleel (|l977h 
and misidentified by iDubner et al.1 (Il99lll as their source 
4, did not have enough informati on to accurately compute a 
spectral index (iMilne et al.lll98^ . Moreover, it was obse rved 
at 80 and 160 MHz with the Culgoora radioheliograph llSled 
Il977h but was no t identified in the Parkes survey at 2.7 GHz 
I Day et'^ll972l ~l. 

Here, we refine the calculation of the spectral indices for 
all five catalogued sources by two methods: T-T plots and a 
linear fit to the logarithmic plot Si, versus v, where v is each 
of the frequencies in which the CABB data were divided, 
and Si, is the flux density of the source. In the first case. 
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-43°30' - 


Right Ascension (J2000) 


Figure 3. Spectral index distribution along Puppis A and surroundings. The filled box of 7x7 pixels at the bottom left corner represents 
the size of the area where T-T plots were computed. All pixels where the linear fits to the T-T plot had a correlation coefficient lower 
than 0.950 were blanked. The bar at the right displays the spectral index. Radio continuum 1.4 GHz contours at 20, 60, 150, 300 and 
500 mjy beam“^ are overlaid. The beam is shown at the top left corner. 


we used the maps filtered with visibilities between 0.5 and 
4 kA (see 1 13.11) . and averaged the spectral indices obtained 
by exchanging both frequencies alternately as abscissa or 
ordinate. To apply the second method, flux densities were 
obtained over the images per individual band (as described 
in the last paragraph of EU by fitting a two-dimensional 
Gaussian to each source and integrating the fitted function, 
allowing for a background level. To ensure the robustness 
of the fits, we made use of two image processing packages: 
MIRIAD and AlPS. Both results were in general agreement, 
and the error bars are consistent with the slight differences 
we measure between the results of these two packages. 

In Table [U we list all five sources and their spectral 
indices as measured with T-T plots (fourth column) and 
with the slope of the Si, versus v curve (fifth column). In all 
cases, the quoted errors correspond to the linear fits, regard¬ 
less of the image sensitivity or other sources of uncertainty, 
like the limits of the integrated regions in the estimation of 
Sv or the background level (we adopted the integrated flux 
density provided by the Gaussian fit over the source in each 
image). The plots and the corresponding fits are displayed 
in Fig. [1 

It is striking that T-T plots show a curvature for sources 
3 and 4, with the highest flux density points departing from 
the fitted line. To illustrate this issue, we try different fits 
to the points and show the results in Fig. [5] The left-hand 
column displays images of sources 3 (top) and 4 (bottom) at 


1.4 GHz constructed with visibilities between 0.5 and 4 kA. 
The right-hand column shows the corresponding T-T plots, 
where the solid lines coincide with the fits shown in Fig. U] 
in which all flux densities above 45 and 35 mJy beam“^ at 

1.4 and 2.5 GHz, respectively for source 3, and above 30 
and 20 mJy beam“^ at 1.4 and 2.5 GHz, respectively, for 
source 4, were taken into account. Limiting the range up to 
100 and 80 mJy beam“^ at 1.4 and 2.5 GHz, respectively, 
for source 3, and up to 80 and 50 mJy beam“^ at 1.4 and 

2.5 GHz, respectively, for source 4, the resulting spectral 
indices are a = —0.5 ± 0.3 for source 3, and a = —1.0 ± 
0.4 for source 4. In Fig. these fits are represented with 
a dashed line. If only the brightest core of each source is 
considered (fluxes above 150 and 120 mJy beam~^ at 1.4 
and 2.5 GHz, respectively, for source 3, and above 80 and 50 
mJy beam“^ at 1.4 and 2.5 GHz, respectively, for source 4), 
the resulting spectral indices are a = —0.14±0.20 for source 
3, and a = —0.55 ± 0.35 for source 4 (dotted lines in Fig. 
EJ. For clarity, the parameters used for all the fits above and 
the corresponding spectral indices derived are summarized 
in Table [5] The curious behaviour depicted by these two 
sources implies that there is a gradient of spectral indices 
from the centre (flatter) to the edges (steeper). 

For completeness, in Table [3] we list all new compact 
sources detected in the field with a flux density above 10 
mJy at 1.3 GHz, estimated on the single 128 MHz band 
image (see last paragraph of E3J. These sources are num- 




[Jy/bean.] ^ Jjy/beam] S, ^ Jjy/beam] ^5 ghz [ Vbeom] 


6 Reynoso and Walsh 










Table 3. New compact sources 



RA (J2000) 
h m s 

Dec. (J2000) 

0 / // 

CO. 

‘^l.SGHz 

mjy 

Spectral 

index 

6 

08 17 58.6 

-43 18 37 

53.5 ± 8.5 


7 

08 18 22.6 

-42 54 12.5 

133 ± 7 

-1.25 ± 0.07 

8 

08 18 35.3 

-43 28 47.0 

180 ± 10 

-1.17 ± 0.08 

9 

08 18 42 

-42 53 03 

20 ± 6 

-0.5 ± 0.4 

10 

08 18 57.9 

-43 31 29.5 

90 ± 7 

-0.78 ± 0.08 

11 

08 19 11 

-42 46 55 

46 ± 6 

-0.75 ± 0.25 

12 

08 19 18.5 

-43 23 55 

10 ± 3.5 


13 

08 19 22.2 

-43 15 09 

76.5 ± 7.5 

-0.70 ± 0.15 

14 

08 19 35 

-43 39 03 

17.5 ± 7.5 


15 

08 19 49 

-43 24 40 

14 ± 4 

-0.1 ± 0.3 

16 

08 19 58 

-43 50 50.3 

49 ± 6 


17 

08 20 25 

-43 47 11 

12 ± 6 


18 

08 20 40 

-42 09 14 

29 ± 4 


19 

08 20 51 

-42 24 00 

18 ± 6 

-0.74 ± 0.30 

20 

08 21 09 

-43 22 35.4 

46 ± 3 


21 

08 21 22 

-42 16 19 

40 ± 6 

-1.12 ± 0.25 

22 

08 21 42 

-42 31 43 

33 ± 8 

-1.0 ± 0.4 

23 

08 21 51 

-43 48 37 

35 ± 6 


24 

08 22 20 

-43 51 24 

18 ± 6 


25 

08 22 38 

-42 30 19 

29 ± 4 

-0.70 ± 0.45 

26 

08 22 56 

-42 26 55 

36 ± 6 

-0.35 ± 0.45 

27 

08 24 27 

-43 00 24 

85 ± 15 


28 

08 24 48 

-43 19 24 

26 ± 5 

-0.53 ± 0.45 

29 

08 24 57 

-43 19 17.3 

37 ± 5 

-0.6 ± 0.3 

30 

08 25 05 

-43 09 21 

12 ± 8 


31 

08 25 06 

-42 42 19 

12 ± 5 


32 

08 25 31 

-42 39 09 

27 ± 6 

-0.8 ± 0.4 

33 

08 25 40 

-43 10 36 

15 ± 4 


34 

08 25 40 

-43 31 29 

18 ± 8 


35 

08 25 53 

-43 24 14 

15 ± 5 


36 

08 26 05 

-43 22 59.4 

100 ± 5 

-1.3 ± 0.5 

37 

08 26 11.7 

-42 39 53.3 

127 ± 8 

-0.9 ± 0.1 

38 

08 26 30 

-43 24 43 

41 ± 5 



^ The flux density at 1.3 GHz is estimated on the single 128 
MHz band image (see last paragraph of EB. 


bered following Table [T] ordered with increasing Right As¬ 
cension. In the fourth column, we give the integrated flux 
density at 1.3 GHz. As in the previously catalogued sources, 
a 2D-Gaussian fit was applied to each source, allowing for 
an offset, and the flux density is computed by integrating 
the fitted distribution. We also include a spectral index esti¬ 
mation based on T-T plots obtained as in Table [T] for those 
sources for which the correlation coefficient of the linear fit 
is better than 0.95. 


4 DISCUSSION 


Figure 4. F 

or each of the five previously known sources, T-T plots 
(left-hand column) and the corresponding linear fits to the flux 
densities S{i') in terms of frequency (right-hand column). The 
numbers at the top right corner of each frame in the right-hand 
column correspond to the number assigned to each source in 
Table [Band in Fig.[B 


The energy distribution of an SNR contains information on 
the mechanisms involved in particle acceleration processes. 
Spectra produced by shock-accelerated radio-synchrotron 
emit ting electrons typically have a spectral index a ~ —0.6 
(e.g. iBietenholz et al.ll2OT]] ). In the test-particle limit, adi¬ 
abatic shocks have a theoretical spectral index a = —0.5, 
assuming a c ompression rati o r = 4 in a monoatomic gas 
(7 = 5/3; e.a'. I RevnoldsI I 2 OTII ) . In radiative shocks, the test- 
particle limit is no longer valid as accelerated particles influ¬ 
ence the shock dynamics, and can increase the compression 
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Figure 5. Top, left-hand panel: compact source 3 at 1.4 GHz. 
The beam is plotted at the bottom left corner. Black contours 
(yellow in the online colour version) at 45, 80, 100, 150 and 200 
mjy beam“^ are overlaid to represent the emission at 1.4 GHz, 
while grey (magenta) contours correspond to the emission levels 
35, 60, 80, 120 and 160 mJy beam“^ at 2.5 GHz. The flux density 
scale, in units of mJy bearn”^, is shown at the right. Top, right- 
hand panel: T-T plot for source 3 with the same fit as in Fig. 
O (solid line) and using different flux density ranges (see text). 
Bottom, left-hand panel: compact source 4 at 1.4 GHz. The beam 
is plotted at the bottom left corner. Grey (magenta) contours at 
30, 50, 80 and 115 mJy beam“^ are overlaid to represent the 
emission at 1.4 GHz, while white (black) contours correspond to 
the emission levels 20, 30, 50 and 75 mJy bearn”^ at 2.5 GHz. 
The flux density scale, in units of mJy bearn”^, is shown at the 
right. Bottom, right-hand panel: T-T plot for source 4 with the 
same fit as in Fig. [4] (solid line) and using different flux density 
ranges (see text). 


ratio and produce flatter spectra. Fig. displays the distri¬ 
bution of spectral indices as obtained from the Ca t alogu e 
of Galactic Supernova Remnants compiled bv iGree nl (|2014h . 
Although the catalogue lists 294 sources, 135 of them have 
variable spectral indices, or an undetermined value (quoted 
value followed by a question mark). In Fig. |6l only the re¬ 
maining 159 SNRs are considered. Spectral indices as ex¬ 
treme as 0.0 or —0.9 are marginally attained, but most SNRs 
have —0.6 a < —0.4. Spectral indices can also be altered 
by contamination with thermal Bremsstralhung emission, 
interaction of the shock front with a denser medium or high 
energetic particles injected by an interior pulsar. 

The most recent determination of the radio spectral in¬ 
dex o f Puppis A has been carried out by ICastelletti et ah] 
ll2006ll who combine a number of flux densities correspond¬ 
ing to frequencies between 19 and 8400 MHz, and obtain 
a = —0.52 ± 0.03. Although in agreement within the er¬ 
ror limits, this value is somewhat flatter than the value in¬ 
ferred here from the T-T plot (Fig. However, we no¬ 
tice that there appears to be a fl attening near 10-100 MHz 
(fig. 4 in ICastelletti et al.l 120091 . where the observational 
points do not follow the fitted single power law model. If 
the two fluxes at the lowest frequencies are excluded, then 
a = —0.56 ± 0.04, in excellent agreement with our result. 

Extended SNRs such as Puppis A provide a unique 
chance to analyse local variations of the spectral index, 
which may reflect local changes in the particle accelera¬ 
tion mechanisms. In general, regio ns of diffuse emission have 
steeper spectra (e.g. Cygnus Loop: [Uvaniker et al.ll20(3i . As 
an extreme case, S147 has a global spectral index of -0.3 but 
in selected diffuse regions it varies between -1.1 and -1.56 
lIXiao et al.ll2004 . The difference between the spectral in¬ 
dices for bright filaments and diffuse regions in S147 is re¬ 
flected in the two branches where the points gather in the 
T-T plots. In Puppis A, in contrast, all points in the T-T 
plot can be fitted by a single line (Fig.[2j, which means that 
this SNR is well represented by a unique global spectral in¬ 
dex. Nevertheless, spatial variations are observed, which will 
be analysed in the following section. 



Spectral Index 


Figure 6. Dis tribution of the spectral indices of catalogued 
Galactic SNRs ([G ree ^ I20l4 . Only 159 out of the 294 catalogued 
SNRs were considered, after discarding variable or unconfirmed 
spectral indices. 


4.1 Spatial spectral index variations 

The first study of the spatial distri bution of spec t ral in - 
dices in Puppis A was performed bv IPubner et ahl l|l99ll l. 
They compare two images, at 327 and 1515 MHz, obtained 
with the Karl G. Jansky Very Large Array (VLA) and inte¬ 
grate the emission in a few boxes over relevant areas. They 
find that the eastern region is considerably steeper than 
the southern, western, northern, and central regions, with 
a =-0.67 ±0.08. 

More recently, ICastelletti et ahl (l2006l i improved the 
analysis by combining new observations at 1425 MHz under¬ 
taken with the VLA with the old 327 MHz data re-processed 
from scratch. To ensure that all scale structures are sampled, 
they integrate single dish observations in their maps. They 
apply two techniques to survey the spectral index distribu¬ 
tion. The first method divides directly the logarithm of the 
ratios between fluxes and frequencies: 
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The second method, called ‘tomography’, consists of scaling 
one of the images by a test spectral index at and subtract¬ 
ing it from the other frequency image. The resulting im¬ 
ages, , will be closer to zero as the test 

spectral index approaches the real one, and features with 
different spectral indices will be highlighted as light (posi¬ 
tive) or dark (negative) as their spectr a are steeper o r flatte r 
than the assumed at- In contrast to iDubner et al.l lll99lll , 
ICastelletti et ah] ll2006tl point out that the bright eastern 
features are somewhat flatter than the rest of the shell. 

Both methods are sensitive to offsets in the background 
emission. In particular, to include the large-scale emission 
at 327 MHz. lCastelletti et al.l (l2006ll use single dish archival 
data at 408 MHz and scale it with a s pectral index a = —0.6 
derived from two im ages, at 86 MHz llMills et al.lll960l 'l and 
at 408 MHz llGreenlflOTlIi . Therefore, the subsequent anal¬ 
ysis can be influenced to some extent by this assumed uni¬ 
form spectral index. Here, we apply a third method which 
is independent of the background emission level and does 
not use archival observations or pre-determined results. The 
method, described in detail in Section 13.11 is called ‘T-T 
plots’, and the spectral index distribution is shown in Fig. 
[3l We first call to attention that this method does not pro¬ 
vide information on the spectral index of the interior diffuse 
emission, since the Mi)-filtering performed to construct the 
maps from a same spatial frequency range in the Fourier 
domain, filters out the underlying large-scale emission. Be¬ 
sides, the size of the region adopted to compute the fit sets 
an additional limit to the scale to which this method is sen¬ 
sitive, since only spectral indices associated with small-scale 
features whose bright ness vary significantly withi n the re¬ 
gion can be measured ll Anderson fc Rudnicklll99)^ . 

To compare our results with previous studies of spec¬ 
tral index distribution, we also calculated T-T plots for the 
regions shown in Fig. [7] For the southern region (bottom 
panel), the T-T plot for the area enclosed by the outermost 
contour, within the selected box, yields a = —0.66 ± 0.20, 
steeper than the value reported bv iDubner et al.l lll99lll . The 
inner peak has a steeper spectrum: a = —0.84:0.4, in agree - 
ment with the tomographic image in ICastelletti et H] (l2006l i 
for at = —0.8, in which most of Puppis A has vanished but 
the southern peak is still visible. T he western ear - like fe ature 
(top, right-hand panel), for which IPubner et al.l (Il99ll l esti¬ 
mated a spectral index of-0.56, could not be reliably fitted: 
we found a = —0.6, but with a very poor correlation coef¬ 
ficient (lower than 0.75). The low emission at this region, 
close to the rms level, must be a reason for the increase 
of the noise of the T-T plot. Besides, as ICastelletti et ahl 
(l2006li noticed, this feature has a combination of different 
structures with variable spectral indices (see also Fig. O- 

The eastern region (top, left-hand panel in Fig. O has 
a global spectral index of a = —0.6 ±0.1 when the emission 
included in the outermost contour for t he whole map iscon- 
sidered. This value is in agreement with iDubner et ^ (1 199 il l 
within the error margins. If the individual boxes are taken 
separately, the results are: asoxi = —0.60 ± 0.45, abox 2 = 
-0.64T0.35, (Tboxs = -0.68T0.25 and abox 4 = -0.55T0.12. 
Curiously, when the T-T plots are limited to the brightest 
regions, the result is different for boxes 3 and 4: for the first 
one, the spectrum steepens (a = —0.8 ± 0.4), while for the 
latter, it gets flatter (a = —0.50 ± 0.35). This means that 




ghgjm^O" Zl“30“ 21“20“ Zl^lO" 21”00“ 

Right Ascension (J2000) 


Figure 7. Radio images of selected regions of Puppis A (top, left- 
hand panel: eastern region; top, right-hand panel: western ear-like 
feature; bottom: southern region) at 1.4 GHz using visibilities 
within the spatial frequency range from 0.5 to 4 kA. The beam 
is shown at the bottom. Representative contours at 1.4 and 2.5 
GHz are overlaid, together with the boxes on which the T-T plots 
were calculated. 


the synchrotron emission and the radio spectral behaviour 
are not necessarily correlated. 

In the spectral index map presented in Fig.|31 T-T plots 
were computed on boxes of ~ 100 x 100 arcsec, which is more 
than twice the beam area. In most cases, regions where the 
spectral indices are more extreme (black and white in Fig.O 
or blue and red in the online colour version) are correlated 
with the larger errors in the spectral index determination 
(Aa ~ 0.6). However, the uncertainties of the flat spectrum 
filaments at the south-eastern region are low enough to en¬ 
sure their relia bility. Our result s do n ot reproduce the stri- 
ation found bv ICastelletti et al.l ll2006li along the shell, with 
fringes perpendicular to the Declination axis. Instead, as 
mentioned above, Fig. [3] shows a pattern of narrow strips at 
the southeastern region parallel to the rim which is not cor¬ 
related with the emission observed in the radio maps (Fig.[T] 
and Mr>-filtered images for T-T plots). To explore the nature 
of these filaments, in Fig. [H] we show an image of Puppis A 
constructed from a more restrictive uv range, with visibili¬ 
ties from 1 to 4 kA, in which extended structures are further 
removed. The eastern region depicts a number of bright fil¬ 
aments whose general appearance resembles the pattern ob¬ 
served in the spectral index distribution, suggesting that the 
variations in a are connected with structures at this scale. 
The fact that t hey do not share a deta iled correlation, as 
no ticed also by Castelletti et al] (120061') in Puppis A and 
bv I Anderson fc RudnickI ( 19931) in the SNRs C39.2-0.3 and 
C41.1-0.3, probably means that the variations in spectral 
index are due to aspects of the hydrodynamic flow unre- 
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Right Ascension (J2000) 

Figure 8. Radio image of Puppis A at 1.4 GHz using visibilities 
within the spatial frequency range from 1 to 4 kA. The beam is 
shown at the bottom left corner. To enhance bright features, the 
intensity is displayed with a logarithmic scale. Contours are as in 

Fig.E] 
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Figure 9. Spectral index versus flux density at 1.4 GHz for Pup¬ 
pis A. Compact sources have been removed. The solid line repre¬ 
sents a linear fit for flux densities above 20 mjy beam“^. 




lated with synchrotron emissivity. Besides, if the spectrum 
of the electron population is intrinsically bent, different de- 
ajrees of compression can re sult in spectral index variations 


grees or compression can re st 
(jAndergg^^^^^ydnicy 199.th 


Anderson fc RudnickI (Il993l 'l mention that in old shell 


fainter or diffuse parts of the remnants (e.g. 

Fuerst & Reichl 

1 19861: lUvaniker et al.ll2004l: lingallinera et al. 

I 2 OI 4 II. while in 


younger SNRs like Puppis A, the brighter regions apparently 
tend to be associated with steeper indices. To see if there 
is evidence of any such trend, in Fig. [5] we plot the spectral 
index against the flux density of Puppis A at 1.4 GHz after 
removing the compact sources. For flux densities above 20 
mJy beam“^, the distribution is fit by an almost horizontal 
line, with a slope of 0.018 ± 0.026. The average spectral in¬ 
dex is —0.57±0.05, with no variations along the flux density 


axis. This result confirms that there is no obvious correlation 
between emissivity and spectral index. 

Flatter spectr a are interpreted as sites o f active par¬ 
ticle acceleration ll Anderson fc RudnickI [l996l i . As a shock 
front encounters gas at a higher density and lower tempera- 
ture, the Mach number increases and the spectr um flattens 
l|Anderson fc R.udnicklll99^ : lLeahv fc Tianl[200Rl '). In Fig. [3l 
it is clear that the filament at the easternmost part of the 
shell has a flatter spectrum than the rest of the rim. This is 
consistent with the presence of an external cloud interacting 
with Puppis A, as suggested by several atomic and molecu- 
lar line stud ies ([Dubner fc Arnailll988l : [R,evnoso et al ] ll99,l 
IParon et al.l 12003 1. In other SNRs, shocked clumps are not 
found to be associated with flat spectral indices (e.g. HB21; 
lLeahvll2006h : this is probably because any interaction with 
a denser gas must be followed by important ionization losses 
for the spectrum to flatten. 

ICastelletti et al.l ll2006l l suggest that the eastern and 
western ‘ears’ observed in radio and X-rays could be the 
termination shocks of two opposite jets ejected by the CCO. 
The spectral index distribution of these features does not 
offer any hint to help conclude if this scenario is true, since 
the eastern ear is clearly flatter than the western one. A 
possible explanation for the difference in a if both features 
share the same origin, would be that the eastern shock has 
encountered a denser ISM and the flatter spectral index is 
due to this interaction. 


4.2 Compact sources 

Previous hi g h-reso l ution radio survey s of Puppis A 
I Milne et al.l 1 19831 : iDubner et al.l Il99li : ICastelletti et al.l 
l2006h have reported the presence of four compact sources 
superimposed on the SNR fac(|3, listed as sources 1 to 4 in 
Table U Source 3 is t he only one w i th an X-ray counterpart 
llDubner et al.|[2013l l. [Milne et al.l (Il983l l estimate spectral 
indices for these four sources comparing their observations 
at 21 cm with the Molongo map at 408 MHz I Green|[l97lll 
and find a ~ —1.0 for all of them. lPubner et al. 1 19911 1 com¬ 
pare flux densities at 327 and 1515 MHz and find spectra 
steeper than —1.0 for all but source 3, for which a = —0.8. 
Base d on the tomographic im age for a test spectral index of 
—0.8. [Castelletti et al.l ll2006h conclude that all four sources 
are steeper than —0.8 and, therefore, extragalactic. Our re¬ 
sults do not confirm these results except for source 2, which 
is clearly extragalactic, together with source 5 which lacks 
previous estimations of the spectral index. However, the re¬ 
sults are more in line if we note that sour ce 3 does not ap¬ 
pear a nymore in the tomographic image of ICastelletti et 51 
(l200fl l with at = —0.8 (hence the spectral index is limited 
between —0.8 < a < —0.6), and source 4 is marginally ob¬ 
served, which means that a ~ —0.8 but not steeper. 

From Table [ll spectral indices measured from the curve 
S(iy) versus u are systematically steeper than those derived 
using T-T plots. For sources 1 and 4, the steeper indices 
would not be in conflict with previous results. However, 


^ A Molonglo contour map at 408 MHz llGree a 11^ shows all 
four sources, but only the brightest one (source 1) is identified 
as an unresolved source with a spectrum steeper than the rest of 
Puppis A. 
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Figure 10. Left: very high resolution (HPBW = 6.8 X 4.6 arcsec) 
image of source 3 at 1.4 GHz. The beam is indicated in the bottom 
left corner. Contours of an image at 1.4 GHz filtered in uv between 
0.2 and 4 kA ( J2]| are overlaid, at 40, 80, 150 and 250 mjy beam“^. 
Right: same as left for source 5, where the plotted contours are 
30, 80, 200, 400 and 500 mJy beam“^. 


spectral indices obtained with the T-T plot method are more 
reliable, since they are free of the following problems in¬ 
volved with flux estimates. In the first place, splitting the 
uv data into 128 MHz bands has the disadvantage that the 
reduced number of visibilities per band increases the noise 
level of the images and makes the removal of sidelobes more 
challenging. And secondly, if the 2D Gaussian or offset fits 
are poor, the integrated flux density may not be accurate. 
The departure of several measured points from the fitted 
lines in the right-hand column of Fig. [4] can be ascribed to 
these two problems. In the case of sources 3 and 4, we have 
discussed in 1)3.21 that a varies from the core to the edges. 
When the core is excluded, the spectral indices derived from 
T-T plots and by fitting a line to SR) versus r are coincident 
(compare Tables [T] and [2|). This means that the contribution 
from the cores is negligible when computing the total fluxes 
of these sources. 

It is possible that the difference between the values of 
Q measured here and in earlier works can be explained be¬ 
cause sources 1-4 are superimposed on the diffuse component 
of Puppis A, and previous flux density determinations are 
probably contaminated with this background emission. It is 
also likely that one or more of these sources are variable. As 
an example, we notice that there is a difference of about a 
factor of 2 between the flux densities of sources 1 and 2 at 21 
cm (500 and 230 mJy, respectively; iMilne et al.lll98 3ll and 
at 15 15 MHz (230 and 160 mJy, respectively; Dubner et al.l 
[1993), whereas according to the quoted spectral indices the 
difference should be less than 10 per cent. In that sense, the 
present work exhibits the unique advantage of being based 
on simultaneous observations. 

Finally, we have made use of the visibilities associated 
with the 6th antenna of the ATCA, which vary in the range 
between 8 and 24 kA, to construct very high resolution im¬ 
ages (HPBW=6.8 X 4.6 arcsec) of the brightest sources. 
These images reveal that sources 3 and 5 are double, as 
shown in Fig. 1101 where the former could be even triple (a 
weaker feature probably real appears at the NW of the two 
brighter components). Such double-lobed structures are typ¬ 
ical of extragalactic sources. The multiple nature of source 3 
could explain the different spectral indices from the edge to 
the centre. However, source 4 does not show any structure 


on small scales, hence the reason of its spectral behaviour 
remains a mystery. Although previously identified as an ex¬ 
tragalactic source, it is more likely a Galactic source that 
may even be associated with Puppis A. 

Among the new compact sources catalogued in Table 
[31 double sources are also: 10, 21, 22, 25, 26 and 38, while 
source 36 is probably triple. Unfortunately, the uv sampling 
obtained with baselines that include only the 6th antenna is 
too scarce to perform a deeper investigation into the prop¬ 
erties of these sources. New observations with the ATGA in 
its 6 km configurations are desirable to uncover the nature 
of these multiple sources. 


4.3 Ejecta ‘bullets’ in Puppis A? 

There is plenty of evidence that core-collapse SN eje ct ‘bul¬ 
lets’ of stellar material (e.g. I Wang fc Ghevali^l200^ . These 
clumps evolve from an initial bow-shock phase to an insta¬ 
bility stage, where Rayleigh-Ta ylor (R-T) fingers a re de¬ 
veloped along the forward face ( Anderson et al.l 1994l 'l. The 
hi gh-velocity optical O -rich knots detec ted in Gassiopeia 
A (Ivan den Berriilll97ll i and G292.0+1.8 ([Ghavamian et ah! 
l2005h are SN ejected clumps in the early bow-shock phase. 
At more evolved stages, R-T fingers stretch and amplify the 
local magnetic field, and the clumps become bright radio¬ 
synchrotron emitters. Compact radio knots are observ ed in 
the young SNR Cassio peia A ([A nderso n et al.l Il994l i and 
in the older Vela SNR ([Duncan et al.lll99fl l. In Gas A, ra¬ 
dio knots have typical spectral indices around —0.75 inde¬ 
pendently of the knot velocity, and at the outer bound- 
aries of the shell, sp ectra are even steeper {a ~ —0.85; 
[Anderson et al.lll994l i. In Vela, several radio k nots are the 
coun t erparts of high-velo city X-ray features (jStrom et al.l 
[199,4 jP uncan et al. Il996li. 

[Gastelletti et al.[ ([2006l i observed the presence of a com¬ 
pact feature (size 45 arcsec) to the east, close to an in- 
tense X-ray regio n known as the ‘brigh t eastern knot’ (BEK; 
[Petre etldl [l982[ : [Hwang et al.[ [2005li . They conclude that 
this source is extragalactic based on a spectrum steeper 
than —0.8 inferred from the tomographic images. In fact, 
this source is hardly observed in their tomographic image 
traced with at = —0.8, which means that the spectral in¬ 
dex is limited between —0.8 < a < —0.6. We analysed our 
images filtered with visibilities between 1 and 4 kA (Fig. 
[8]) and could locate this source at RA(2000)=8^24“'13'!9, 
Dec.(2000)=—42°59'55". The source is resolved, and its size 
is about 86 x 42 arcsec after deconvolving the fitted Gaus¬ 
sian by the beam. A T-T plot yields a = —0.8 ± 1.1. Thus, 
not only is the proposed extragalactic nature of this source 
not unequivocal, but the spectral index is compatible with 
a radio knot, like those found in Gas A. 

Remar kably, this compact radio feature is located next 
to the BEK. lilwang et ^ (|2005[ ') interpreted the X-ray emis¬ 
sion from the BEK as coming from a shocked interstellar 
cloud based on the ‘voided sphere’ morphology similar to 
that produced in experimental simulations of interaction be¬ 
tween shock waves and dense clouds. However, a molecu¬ 
lar study failed to detect any gas concentration associated 
with the BEK ([Paron et al. [2008[ l. It should not be dis- 
carded that both the B EK and the radio source reported by 
[Castelletti et al.) ([200(4 are evolved ejecta bullets like those 
detected in X-rays and radio wavelengths in other SNRs. 
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Table 4. Radio knot candidates® 


RA (J2000) 
h m s 

Dec. (J2000) 

o / // 

Spectral 

index 

08 21 38.3 

-42 38 22 

-0.8 ± 0.5 

08 21 51.4 

-42 39 51 

-0.7 ± 0.7 

08 22 12.5 

-43 19 36 

-0.6 ± 1.1 

08 22 50.4 

-42 39 09 

-0.8 ± 0.9 

08 23 18.7 

-42 52 34 

-0.8 ± 1.1 

08 24 12.4 

-43 09 58 

-0.6 ± 1.0 

08 24 13.9 

-42 59 55 

-0.8 ± 1.1 

08 21 42.2 

-42 41 09 

-1.3 ± 0.6 

08 22 59.4 

-42 36 06 

-1.1 ± 0.3 

08 23 07.7 

-42 42 24 

-1.3 ± 0.5 

08 23 14.1 

-42 41 56 

-1.2 ± 1.2 

08 23 47.3 

-42 51 33 

-1.2 ± 1.6 

08 21 56.0 

-43 20 51 

-0.2 ± 0.7 

08 23 56.37 

-43 12 09 

-0.35 ± 0.60 


“ Radio knot candidates are listed in the upper band. For com¬ 
pleteness, sources probably extragalactic found in the search are 
included in the central band, and flat spectrum sources, in the 
lower band. 

Moreover, the coexistence of ejecta clumps in their earlier 
stages (as optical O-rich knots) and in their more evolved 
phase (as radio knots) in Cas A suggests that radio knots 
can also be found in Puppis A, since both remnants share the 
characteristic system of O-rich high-velocity clumps ejected 
in the SN event. 

We have searched for radio knots along the shell and 
found six more candidates, based on their spectral indices, 
which are displayed in Table U) We also list five compact 
sources with spectra steeper than -1, hence probably of ex¬ 
tragalactic origin (although the error margins make their 
classihcation not conclusive), and two compact sources with 
a flat spectrum. Higher resolution and sensitivity images of 
these radio knot candidates will help to measure changes 
in intensity and position, and to unveil if they are individ¬ 
ual clumps or are composed by a cluster of knots. Proper 
motion studies are crucial to conclude if these features are 
high-velocity ejecta. 


5 CONCLUSIONS 

Using the ATCA, we have imaged the prominent SNR Pup¬ 
pis A in continuum emission between 1 and 3 GHz. The 
observations included full 12 hr synthesis imaging in two ar¬ 
rays - EW 352 and 750 A. Our images are composed of 24 
mosaic pointings and cover an area approximately 2x2 deg^. 
The resultant data are of high sensitivity (1.5 mjy beam“^) 
and high resolution (~ 80 x 50 arcsec). 

We have used T-T plots to calculate the spectral in¬ 
dex across Puppis A and derive a value of —0.563 ± 0.013. 
Although local variations in the spectral index are noted, 
this global spectral index appears to well represent most 
of the SNR. Contrary to previous suggestions, we End no 
correlation between spectral indices and radio brightness. 
We identify filamentary structure, particularly in the east¬ 
ern edge of the SNR, which resembles the pattern followed 
by the spectral index variations. This pattern was not de¬ 


tected in previous studies. We suggest that the outermost 
filament may be the result of interaction of the SNR with an 
external cloud. There is no evidence, based on the spectral 
index, of the eastern and western ‘ears’ being termination 
shocks of two opposite jets from the CCO. We also found 
a number of compact features that could be radio knots re¬ 
sulting from the evolution of clumpy SN ejecta. A dynamical 
study should clarify this hypothesis. 

We investigate the nature of a number of compact 
sources that are found projected both within and outside of 
the SNR. Some of these compact sources do show strongly 
negative spectral indices, which are indicative of extragalac¬ 
tic sources. It is not clear that sources 1, 3, and 4 are ex¬ 
tragalactic, as claimed in previous works. With the highest 
spatial resolution available to us, we are able to identify 
double-lobed structures in sources 3 and 5, confirming the 
extragalactic nature of the latter and suggesting the same 
for the former in spite of its flat spectral index. However, 
source 4, which was previously identified as extragalactic, 
shows a flatter spectral index with a radial gradient (like 
source 3) and does not show any structure on small scales. 
We therefore expect that source 4 is more likely a Galactic 
source that may well be associated with the SNR. 


ACKNOWLEDGEMENTS 

We thank the referee (Luke Bozzetto) for helpful and con¬ 
structive criticism, which has improved the quality of this 
paper. We appreciate useful discussions with Laura Richter, 
Timothy Shimwell and Mark Wieringa for solving problems 
during data reduction. This research was partially funded by 
CONIGET grants PIP 114-200801-00428 and 112-201207- 
00226. The Australia Telescope Compact Array is part of 
the Australia Telescope National Facility which is funded 
by the Commonwealth of Australia for operation as a Na¬ 
tional Facility managed by CSIRO. EMR is member of the 
Carrera del Investigador Cientfflco of CONIGET, Argentina. 


REFERENCES 

Anderson M. C., Jones T. W., Rudnick L., Tregillis 1. L., 
Kang H., 1994, ApJ, 421, L31 

Anderson M. C., Rudnick L., 1993, ApJ, 408, 514 

Anderson M. C., Rudnick L., 1996, ApJ, 456, 234 

Arendt R. G., Dwek E., Blair W. P., Ghavamian P., Hwang 
U., Long K. S., Petre R., Rho J., Winkler P. F., 2010, ApJ, 
725, 585 

Arendt R. G., Dwek E., Petre R., Dickel J. R., Roger R. S., 
Milne D. K., Kesteven M. J., 1990, ApJ, 350, 266 

Becker W., Prinz T., Winkler P. F., Petre R., 2012, ApJ, 
755, 141 

Bietenholz M. F., Kassim N. E., Weiler K. W., 2001, ApJ, 
560, 772 

Castelletti G., Dubner G., Golap K., Goss W. M., 2006, 
A&A, 459, 535 

Costain C. H., 1960, MNRAS, 120, 248 

Day G. A., Caswell J. L., Cooke D. J., 1972, Aust. J. Phys. 
Astrophys. Suppl., 25, 1 

Dickel J. R., Milne D. K., 1976, Aust. J. Phys., 29, 435 








12 Reynoso and Walsh 


Doherty M., Johnston S., Green A. J., Roberts M. S. E., 
Romani R. W., Gaensler B. M., Crawford F., 2003, MN- 
RAS, 339, 1048 

Dubner G., Loiseau N., Rodriguez-Pascual P., Smith 
M. J. S., Giacani E., Castelletti G., 2013, A&A, 555, A9 
Dubner G. M., Arnal E. M., 1988, A&AS, 75, 363 
Dubner G. M., Braun R., Winkler P. F., Goss W. M., 1991, 
AJ, 101, 1466 

Duncan A. R., Stewart R. T., Haynes R. F., Jones K. L., 
1996, MNRAS, 280, 252 
Fuerst E., Reich W., 1986, A&A, 163, 185 
Gaensler B. M., Brazier K. T. S., Manchester R. N., John¬ 
ston S., Green A. J., 1999, MNRAS, 305, 724 
Ghavamian P., Hughes J. P., Williams T. B., 2005, ApJ, 
635, 365 

Green A. J., 1971, Aust. J.Phys., 24, 773 
Green D. A., 1990, AJ, 100, 1927 
Green D. A., 2014, Bull. Astron. Soc. India, 42, 47 
H. E. S. S. Collaboration : Abramowski A., Aharonian F., 
Ait Benkhali F., Akhperjanian A. G., Angiiner E. O., 
Backes M., Balenderan S., Balzer A., et al. 2015, A&A, 
575, A81 

Hewitt J. W., Grondin M.-H., Lemoine-Goumard M., Re¬ 
poseur T., Ballet J., Tanaka T., 2012, ApJ, 759, 89 
Hui C. Y., Becker W., 2006, A&A, 457, L33 
Hwang U., Flanagan K. A., Petre R., 2005, ApJ, 635, 355 
Hwang U., Petre R., Flanagan K. A., 2008, ApJ, 676, 378 
Ingallinera A., Trigilio C., Umana G., Leto P., Agliozzo C., 
Buemi C., 2014, MNRAS, 445, 4507 
Katsuda S., Hwang U., Petre R., Park S., Mori K., Tsunemi 
H., 2010, ApJ, 714, 1725 

Katsuda S., Mori K., Tsunemi H., Park S., Hwang U., Bur¬ 
rows D. N., Hughes J. P., Slane P. O., 2008, ApJ, 678, 297 
Kundu M. R., 1970, ApJ, 162, 17 
Leahy D., Tian W., 2005, A&A, 440, 929 
Leahy D. A., 2006, ApJ, 647, 1125 

McConnell D., Sadler E. M., Murphy T., Ekers R. D., 2012, 
MNRAS, 422, 1527 

Mills B. Y., Slee O. B., Hill E. R., 1960, Aust. J.Phys., 13, 
676 

Milne D. K., 1971, Australian Journal of Physics, 24, 429 
Milne D. K., 1972, Aust. J.Phys., 25, 307 
Milne D. K., Goss W. M., Danziger 1. J., 1983, MNRAS, 
204, 237 

Milne D. K., Hill E. R., 1969, Aust. J.Phys., 22, 211 
Milne D. K., Stewart R. T., Haynes R. F., 1993, MNRAS, 
261, 366 

Paron S., Dubner G., Reynoso E., Rubio M., 2008, A&A, 
480, 439 

Petre R., Becker C. M., Winkler P. F., 1996, ApJ, 465, L43 
Petre R., Kriss G. A., Winkler P. F., Canizares C. R., 1982, 
ApJ, 258, 22 

Reich P., Reich W., 1988, A&AS, 74, 7 
Reynolds S. P., 2001, Space Sci. Rev., 99, 177 
Reynoso E. M., Dubner G. M., Goss W. M., Arnal E. M., 
1995, AJ, 110, 318 

Reynoso E. M., Green A. J., 2007, MNRAS, 375, 92 
Sault R. J., Teuben P. J., Wright M. C. H., 1995, in 
R. A. Shaw, H. E. Payne, & J. J. E. Hayes ed.. Astro¬ 
nomical Data Analysis Software and Systems IV Vol. 77 
of ASP Conf. Ser., A Retrospective View of MIRIAD. pp 
433—f 


Slee O. B., 1977, Aust. J. Phys. Astrophys. SuppL, 43, 1 
Strom R., Johnston H. M., Verbunt F., Aschenbach B., 
1995, Nature, 373, 590 

Turtle A. J., Pugh J. F., Kenderdine S., Pauliny-Toth 
1. 1. K., 1962, MNRAS, 124, 297 
Uyamker B., Reich W., Yar A., Fiirst E., 2004, A&A, 426, 
909 

van den Bergh S., 1971, ApJ, 165, 457 
Wang C.-Y., Chevalier R. A., 2002, ApJ, 574, 155 
Watson M. G., Schroder A. C., Fyfe D. e. a., 2009, A&A, 
493, 339 

Winkler P. F., Kirshner R. R, 1985, ApJ, 299, 981 
Winkler P. F., Petre R., 2007, ApJ, 670, 635 
Winkler P. F., Tuttle J. H., Kirshner R. P., Irwin M. J., 
1988, in Roger R. S., Landecker T. L., eds, lAU Colloq. 
101: Supernova Remnants and the Interstellar Medium 
Kinematics of Oxygen-Rich Filaments in Puppis A. p. 65 
Woermann B., Gaylard M. J., Otrupcek R., 2000, MNRAS, 
317, 421 

Xiao L., Fiirst E., Reich W., Han J. L., 2008, A&A, 482, 
783 

Zavlin V. E., Triimper J., Pavlov G. G., 1999, ApJ, 525, 
959 

This paper has been typeset from a iQijX/ DTf^X file prepared 
by the author. 


